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We review what is different and what is similar in a color superconductor as com- 
pared to an ordinary BCS superconductor. The parametric dependence of the 
zero-temperature gap, </>o, on the coupling constant differs in QCD from that in 
BCS theory. On the other hand, the transition temperature to the superconduct- 
ing phase, T c , is related to the zero-temperature gap in the same way in QCD as 
in BCS theory, T C /0 O 0.567. 



1 Cold, dense fermionic matter 

Consider a degenerate, non-interacting fermionic system, where all momentum 
states up to the Fermi surface are occupied. The fcrmion occupation number 
is a step function, no(eo) = 9(— eo), where eo = u> — \i is the energy of the 
fermions with respect to the Fermi energy, fx, and u> is their kinetic energy. 
For non-relativistic particles with 3-momentum k and mass m, u> = k 2 /(2?n), 
while for massless (ultrarelativistic) particles, u> = k = |k|. (The units are 
h = ks = c = 1.) The excitation spectrum of the fermions consists of a 
particle branch, 6q = eo, and a hole branch, e{j = — eo- The formation of 
particle-hole excitations at the Fermi surface costs no energy, Cp + e{j = 0. 

Now switch on an interaction, for the sake of simplicity a point-like four- 
fcrmion interaction with interaction strength G 2 . Let the sign of G 2 be defined 
such that G 2 > in the case of an attractive interaction, and G 2 < 0, if 
the interaction is repulsive. Due to Pauli's exclusion principle and energy 
conservation, scattering of fermions can occur exclusively at the Fermi surface. 
In other words, physical scattering processes, are possible only for eo — >■ 0. The 
amplitude for fermion-fermion scattering isa 

r(£0) ~ 1 - GMn( M /eo) ' (1) 
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where we suppress all constant factors that are irrelevant for the qualitative 
arguments presented here. As eo — > 0, in the case of an attractive interaction 
the scattering amplitude develops a singularity at an energy scale 

e* ~ /i exp(-l/G 2 ) . (2) 

For repulsive interactions, no such singularity occurs. Obviously, the sin- 
gularity (^) occurs even when the attractive interaction is arbitrarily weak, 
G 2 — * + , all that changes is the energy scale £q of the singularity. 

This singularity is the famous Cooper instability El. It is cured by the for- 
mation of Cooper pairs which, as bosons, condense in the true, energetically 
favored ground state of the system. Macroscopically, this leads to the phe- 
nomenon of superfluidity, or, if the fermions carry charge, superconductivity. 

The excitation spectrum of the system also changes. A gap 0o forms at 
the Fermi surface, separating the branch for quasiparticle excitations, e p = — e, 
e = yig + </>q , from the one for quasihole excitations, e h = +e. Now exciting 
a quasiparticle-quasihole pair costs at least an energy e p + e h > 2 0o- The 
quasiparticle occupation number, n(eo), is "smeared" around the (original) 
Fermi surface, n(eo) = (e — eo)/(2e). n 

To compute the gap, one has to solve a gap equation El. In the case of a 
point-like four-fermion interaction, this equation takes the form 

4> = G* f ^<j> . (3) 

(Again, irrelevant constant factors are omitted.) Besides the trivial (ener- 
getically disfavored) solution (f>o = 0, this equation has always a non-trivial 
(energetically favored) solution, ^ 0. Since the gap is a constant for point- 
like four-fermion interactions, we can divide both sides of Eq. (0) by <pQ. The 
remaining integral can be solved exactly, with the result 

0o~2 A iexp(-l/G 2 ) . (4) 

Apparently, (f>o > e^, cf. Eq. (||). In other words, since there are no quasiparticle 
states with energy e p > — cf>o (or quasihole states with e h < 0o), the gap has 
just the right order of magnitude to prevent the scattering amplitude (Q) from 
developing a singularity. 

The assumption of a point-like four-fermion interaction can be relaxed. 
Assume that the interaction is mediated by a scalar boson of mass M. For 
the sake of definiteness, assume that the boson mass is generated by many- 
body interactions at nonzero density, M ~ gfj,. The boson-fermion coupling is 
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denoted by g, the boson propagator is A(P) = 1/(M 2 — P 2 1 P = P^ = (po,p), 
P 2 = p 2 — p 2 . In this case, the gap equation (pi) becomes El 



4>a{k) = g 2 J -^r-(/>o(q) | ha 



M 2 + {k + qf 



M 2 + {k - q) 2 



(5) 



Here, q = eg + /i. The frequency dependence of the boson propagator and, 
consequently, that of the gap function has been neglected, based on the argu- 
ment that in weak coupling, g -C 1, Po ~ e ~ <fio ~ A* ex P(~ 1/.9 2 ) "C M ~ 17/1. 
The logarithm arises from the integration over the angle between the boson 
3-momentum p = k— n and the 3-momentum of the fermion in the condensate, 
q, for details see Ref. □. Note that this factor enhances contributions from the 
region of momenta q ~ k (collinear enhancement). 

In the case of a boson-mediated interaction, the gap is no longer constant, 
but a function of momentum. Consequently, the gap equation (SI) is no longer 
a simple fix-point equation for 0o, but an integral equation which has to be 
solved numerically. However, to estimate the order of magnitude of the gap at 
the Fermi surface, k = \i (for massless fermions), one may make the following 
approximation. First note that, due to the factor 1/e, the integrand peaks 
at the Fermi surface. It is then sufficient to approximate the slowly varying 
logarithm and the gap function with their values for q = k = fi. This leads to 
the estimate 

1 



2 \x exp 



g 2 ln(l + 2^ 2 /M 2 ) 



(G) 



for the value of the gap function at the Fermi surface, </>o = <po(fi), which should 
be compared with Eq. (H|). All that changed is that the coupling constant is 
effectively increased by the logarithm originating from collinear enhancement. 
For M ~ gn, the logarithm becomes ~ ln(l/</) in weak coupling. 



2 Cold, dense quark matter 

The density in cold quark matter increases ~ /i 3 . Asymptotic freedom then 
implies that single-gluon exchange becomes the dominant interaction between 
quarks. Single-gluon exchange is attractive in the color-antitriplet channel, 
and. therefore leads to color superconductivity in cold, dense quark matter 
3 C'EI. Recently, considerable activity was generated 6 ~ 39 by the work of Refs. 
3d, which suggested that the zero-temperature color-superconducting gap (f>o 
could be as large as 100 MeV. This order of magnitude is quite surprising, 
because earlier work by Bailin and Love El estimated the gap to be 4>q ~ 1 
MeV. While gaps of order 100 MeV could also be relevant for the physics 
of nuclear collisions (see below), gaps of about 1 MeV allow at most for the 
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possibility that neutron star cores, if consisting of quark matter, could be color 
superconductors. __ 

The authors of Refs.Q'EI based their arguments on a simple model where 
quarks interact via a point-like four-fermion interaction, giving rise to a gap 
equation of the type ([|). The coupling strength G 2 was adjusted such that the 
model reproduced the order of magnitude of the chiral transition at nonzero 
temperature and zero quark chemical potential, T x ~ 150 MeV. The earlier 
work of Bailin and LoveB was already more sophisticated in the sense that 
they used one-gluon exchange, employing gluon propagators with electric and 
magnetic screening masses. This gives rise to a gap equation of the form (|J), 
with g being the strong coupling constant. 

Unfortunately, both approaches fail to capture an essential property of 
single-gluon exchange: at zero temperature, due to the absence of magnetic 
screening, magnetic interactions are truly long-range. Surprisingly, this fact 
was already known to Barroista, but apparently never made it into the published 
literature. 

Long-range magnetic interactions have the important consequence that 
one can no longer neglect the frequency dependence of the boson propagator, 
as done in the derivation of For massless boson exchange, the gap equation 
assumes the (approximate) form 

0o(fco)-.9 2 /°^0o( e )iln(^^) • (7) 

Again, as in (^) there is an additional logarithm, representing collinear en- 
hancement. In QCD it arises from the exchange of ultrasoft, magnetic gluons. 
While the collinear enhancement in Eq. (Q) is cut off by the mass of the scalar 
boson, M, here it is cut off by the energy of the magnetic gluon, po = eifco- For 
e ~ fco ~ 0o i Po is of the order of <ji>o, too, while in weak coupling, M ~ gy, 3> 4>q. 
For magnetic gluon exchange in QCD, the contribution of the collinear region, 
q ~ k, to the gap integral is therefore much larger than in the case of massive 
boson exchange. 

To estimate the effect of the logarithm on the solution of the gap equation, 
let us neglect the energy dependence of the gap function in the integrand and 
consider its value at the Fermi surface, k = fi, ko = 4>o- One may also make 
the approximations ln[/i 2 /(e 2 — 0q)] ~ 2 ln(/i/e). Then, the integral is again 
exactly solvable, with the (order of magnitude) result 

4>o^2 M e X p(-l/ 5 ) . (8) 

Due to the explicit e dependence of the logarithm in (]?]), the power of g in 
the exponent is reduced as compared to the BCS result (0). The case of 
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massive scalar boson exchange, Eq. (Q), interpolates between these two cases, 
as g 2 <C g 2 ln(l/g) <C g for g <C 1. This reflects the fact mentioned earlier 
that, while the collinear singularity q = k in (||) is cut off by the mass of the 
scalar boson, M ~ gfi, the singularity in (0) is cut off by the gluon energy, 
Po ~ (f> 0i which is, in weak coupling, much smaller than M. In the literature, 
the parametric dependence. nn g of the solution (||) to the gap equation (Q) was 
first discussed in Refs. oElo. 

A^nf Ipday, various refinements of the solution (||) have been discussed 
The value of the gap function at the Fermi surface, k = (i, 
ko = 0o, is 



= b n g 5 exp 



-I) [1 



0(g)] , c=-= 



3_vr2 
V2 



(9) 



Furthermore, the gap function has a non-trivial (on-shell) energy dependence, 



4>o (e) = 0o sin 



6 // 

W) 



(10) 



The constant c was first computed by Son0. To obtain the correct numerical 
value for c, one has to account for-the modifications of the gluon propagator 
in the presence of a dense medium c2l. Then, what dominates the gap equation 
(Q) and determines c is the contribution from nearly static, Landau-damped 
magnetic gluons. Furthermore, Son showed that in computing c, it is essential 
to retain the energy dependence (|l^) of the gap function. To the level of 
accuracy considered by SoncJ, b — 1 and /(e) = e. The prefactor g~ b arises 
from subleading contributions of static electric and non-static magnetic gluons. 

The constant b collects constant factors inJdiese subleading contributkms. 
It was first computed by Schafer and WilczekEJ and the present authors t3, 



5/2 
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(11) 



rafer and WilczekEl obtained b' = 1/2 
that actually b' = 1 and /(e) = e — £q- 



b = 512 7T 4 

with an undetermined constant 6'. f 
and, as before, /(e) = e. We showed 

The additional factor 2 is the same that occurs in Eqs. (|J), (||), and (|8|). It arises 
from the measure of integration in the gap equation, deo/e = — dln(e — eo), 
and — ln(e — eo)! ^ ~ ln(2ui0o)- This also explains the modification of /(e). 
In addition, we pointed outEj that it is important to compute the gap function 
on the correct quasiparticle mass shell, and we considered for the first time the 
case of non-zero temperature. As will be discussed in more detail below, the 
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a„: 3-loop, N,=3, A=364 MeV 

200 | . . . . 1 . . . . r- 



N f =2 




[i [MeV] 

Figure 1: The value of the zero-temperature gap at the Fermi surface, </>o, as a function of 
the quark chemical potential /i for 2 quark flavors (full line), 3 quark flavors (dotted line), 
and 3 quark flavors including effects from the quark wavefunction renormalization in the gap 
equation (dashed line). 



temperature T c where the color-superconducting condensate melts is related to 
the zero-temperature gap in the same way as in BCS theory, T c /(f>o = e 7 /tt ~ 
0.567, where 7 ~ 0.577 is the Euler-Mascheroni constant. 

Brown, Liu, and Ren E3 computed T c from -the quark-quark scattering 
amplitude. Using the result T c /4>o — e 1 j-K of Ref.Ej, they obtain 

b' = exp ( ~ 0.176 , (12) 

where the correction to the previous estimate b' = 1 arises from a finite, \i 
dependent contribution to the wavefunction renormalization for quarks in a 
dense medium. The authors of Ref. E3 also assert that there are no further 
corrections to b' at this order in [f correct, this is a remarkable result, 
because from previous calculations EjIIj it appeared that computing b' exactly 
to leading order in g would be a formidable task. 

Unlike calculations of the free energy or the Debye mass, wiiere the per- 
turbative expansion in powers of g 2 appears to be well-behaved ta even when 
extrapolated down to moderate values of (A, this result for the wavefunction 
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renormalization (which is equivalent to non-Fermi liquid behavior) indicates 
that perturbation theory is not such a good approximation, at least for quarks 
near the Fermi surface. The confirmation of the results of Ref.EZI is clearly an 
outstanding problem for the field. 

Figure [l] shows the value of the zero-temperature gap <po at the Fermi 
surface according to Eq. @ as a function of [i for Nf — 2 and 3 massless 
quark flavors with b' = 1 (full and dotted lines), and for Nf = 3 flavors with 
b' = cxp[— (it 2 + 4)/8] (dashed line). The running of the coupling g(/i) with 
the chemical potential /i was computed from the 3-loop QCD (3 functioned, 
however, not for 6 but only for 3 flavors of massless quarks. Therefore, the 
QCD scale A = 364 MeV is chosen somewhat larger than the standard value, 
to give the value a s (/i = 2GeV) ~ 0.309. Although an extrapolation of the 
weak-coupling result (||) to large g (small /x) appears audacious, it is interesting 
to note that the maximum value of the gap for b' = 1 is of the.-prder of 100 
MeV, quite in agreement with the earlier estimates of Refs. Qu. However, 
taking into account the quark wavefunction renormalizationEZI reduces the gap 
to values of a few MeV. These values are of the order of typical superfiuid 
gaps in ordinary hadronic matter. This lends credibility_to the conjecture that 
quark and hadronic matter are continuously connected^, although symmetry 
arguments L3ii3 suggest that, at zero temperature, there is a first order phase 
transition between these two phases of nuclear matter. 



3 Not so cold, dense quark matter 

To understand how the color-superconducting gap changes with temperature, 
it is instructive to first consider the simpler BCS case. At nonzero temperature 
T, the gap equation (||) becomes 

= G 2 f° deo tanh 

Here, </> is the value of the gap at temperature T, <p = <f>(T), and as before, tf>o = 
0(0) denotes the zero-temperature gap in the following. The <f> dependence of 
the quasiparticle excitation energy e((f>) = y^eg + 4> 2 has been made explicit, 
to distinguish e(</>) from e = e(</>o) used previously. 

Again, (j) is constant and can be divided out on both sides of Eq. ([l3]) , with 
the result 

The effect of the tanh is to reduce the value of the integrand, such that <fi m 
the factor l/t{4>) has to decrease in order to balance the 1 on the left-hand 



2T 



(13) 



2T 



(14) 
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side. At some critical temperature T c , this balance can no longer be achieved 
and 4> = is the only solution of Eq. (p"3|). T c is the temperature where the 
superconducting condensate melts. Physically, the random thermal energy of 
the fermions exceeds their binding energy in a Cooper pair. Thus, T c must be 
of the same order as 0o- 

There is an easy way to comniite the change of </> with T, which was to our 
knowledge first suggested in Ref. 15. Note that tanh[e(</>)/2T] ~ 1 far from the 
Fermi surface, e(</>) > | eo | 3> <Po ~ T. A nonzero temperature influences the 
integrand only in the region close to the Fermi surface, |eo| < 4>o ~ T. Let us 
therefore divide the range of integration into two parts, > eo > —ntfio, and 
—k4>q > 6q > —fi, where k> 1. Then, the tanh need only be kept in the first 
region, and the integral over the second region can be performed similarly as 
at zero temperature, 



1 ~ G 



de 

<4>) 



tanh 



<4>) 



2T 



G In 



K(/)0 



(15) 



Using the solution (Q) for the zero-temperature gap <po, the second term be- 
comes 1 — G 2 ln(2/t). Apparently, the 1 on the left-hand side is almost com- 
pletely saturated by this second term. Cancelling the 1 and writing ln(2 n) = 
f_ K ± deo/e(^o)j one obtains the condition 



G 



de 



K 4>0 



<<p) 



tanh 



2T 



e(0o) 



= 



(16) 



The dependence on k is spurious: one might as well send n — > ooEl, because 
the integrand vanishes when k ^> 1. Equation ( |l^ ) determines <j>(T)/<f>o as a 
function of T. In particular, at T c , where <f> = 0, one derives the well-known 
result 

4^ = — ^ 0.567 (17) 



mentioned above. 

In QCD, it turns outlij that the effect of temperature on the gap equation 
is essentially identical to that in BCS theory: the integrand in (f|) is multiplied 
with tanh[e(0)/2T]. For the same reasons as in the BCS case, this factor is 
negligible far away from the Fermi surface. One may again divide the range 
of integration into two parts, and neglect the effects of temperature in the one 
far from the Fermi surface. The integral over this region can be computed as 
for T = 0. Quite similarly to the treatment in the BCS case, it is found to 
saturate the left-hand side of the gap equation up to corrections of order O(g). 
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Therefore, the integral over the region close to the Fermi surface must 
also be of order O(g), in order to cancel these corrections. To see this, it is 
permissible to compute this integral to leading order in g. One then derives the 
samp_£ondition (|l6| ) as in the BCS case, except that G 2 is replaced by g, see 
Ref.ll3 for details. Consequently, to leading order in g, the T dependence of the 
gap at the Fermi surface, normalized to the zero-temperature gap, </>(T)/</>o, 
is the same as in BCS theory. In particular, the ratio T c /0o is again given 
by Eq. (17). In retrospect, this is not surprising: the prefactor b of the zero- 
temperature gap, Eq. (pi]), was seen to be determined by subleading terms 
in the gap equation [terms of order 0(g) relative to the leading terms due to 
Landau-damped magnetic gluons]. As explained above, temperature affects 
the gap equation at the same subleading order. 

An immediate consequence is that when multiplying the ordinate of Fig. 
1 with 0.567, one obtains the location of the phase transition to the color- 
superconducting phase in the T — fi phase diagram of nuclear matter. In the 
2-flavor case without wavefunction corrections, the transition temperature is 
of order ~ 100 MeV. The color-superconducting phase could then be accessible 
in heavy-ion collisions at BNL-AGS or GSI-SIS energies, which explore the 
range of moderate temperatures and high (net) baryon density in the nuclear 
matter phase diagram. However, in the 3-flavor case, including the effects of 
the quark wavefunction renormalization, the transition temperature is at most 
~ 6 MeV. For such small temperatures, color superconductivity occurs at best 
in neutron star cores, if they consist of quark matter. 
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